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We investigate a two-species Fermi gas with a large mass ratio interacting by an interspecies 
short-range interaction. Using the Born-Oppenheimer approximation, we determine the interaction 
energy of two heavy fermions immersed in the Fermi sea of hght fermions as a function of the s- 
wave scattering length. In the BCS limit, we recover the perturbative calculation of the effective 
interaction between heavy fermions. The p-wave projection of the effective interaction is attractive 
in the BCS limit while it turns out to be repulsive near the unitarity limit. We find that the p- 
wave attraction reaches its maximum between the BCS and unitarity limits, where the maximal 
p-wave pairing of heavy minority fermions is expected. We also investigate the case where the heavy 
fermions are confined in two dimensions and the p-wave attraction between them is found to be 
stronger than that in three dimensions. 
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I. INTRODUCTION 

Experiments using ultracold atomic gases have 
achieved great success in realizing a new type of fermionic 
superfluid. By arbitrarily varying the strength of interac- 
tion via the Feshbach resonance, the crossover from the 
BCS superfluid of fermionic atoms to the Bose-Einstein 
condensate of tightly bound molecules has been observed 
and extensively studied [1, 2]. Nowadays a portion of 
interests of the cold atom community is shifting to the 
two-species Fermi gas with unequal densities [3-7] and 
with unequal masses [8-10]. Elucidating the phase dia- 
gram of such an asymmetric Fermi gas is an important 
and challenging problem. Such an asymmetric system of 
fermions will be also interesting as a prototype of high 
density quark matter in the core of neutron stars, where 
the density and mass imbalances exist among different 
quark flavors [11]. 

So far several quantum phases have been proposed as a 
ground state of the density-imbalanced Fermi gas. One of 
such phases is the p-wave superfluid phase by the pairing 
between the same species of fermions [12]. In Ref. [12], 
the attraction between the same species of fermions was 
found to be induced by the interaction with the other 
species of fermions based on the perturbative calculation 
in the weak-coupling BCS limit. Although the resulting 
pairing gap is exponentially suppressed in the BCS limit, 
it may be possible that the pairing gap becomes large 
enough near the strong-coupling unitarity limit. Thus 
an important question is how large the pairing gap can 
be near the unitarity limit. 

However, difficulties for theoretical treatments arise 
away from the weak-coupling BCS limit because we do 
not have a controlled tool to analyze the strongly in- 
teracting many-body system. Monte Carlo simulations 
also have limitations in asymmetric systems due to the 
fermion sign problem. We note that in one dimen- 
sion the Monte Carlo simulations do not suffer from 
the sign problem and have been employed to study 
the Fulde-Ferrel-Larkin-Ovchinikov phase in the density- 



imbalanced Fermi gas [13-15]. 

In order to obtain further insight into the density- 
imbalanced Fermi gas and the p-wave pairing therein, 
we investigate the Fermi gas with unequal masses be- 
tween two different atomic species. In the limit of large 
mass ratio, we can determine the effective interaction 
among heavy fermions induced by the interaction with 
light fermions using the Born-Oppenheimer approxima- 
tion. Because we do not need to rely on the weak- 
coupling approximation, we can go beyond the pertur- 
bative BCS regime to the strongly interacting unitary 
regime in a controlled way. 

We will show that the p-wave projection of the effec- 
tive interaction between two heavy fermions is attractive 
in the BCS limit being consistent with the perturbative 
prediction, while it turns out to be repulsive near the 
unitarity limit. Thus the p-wave attraction has a max- 
imum between the BCS and unitarity limits, where the 
maximal p-wave pairing of heavy minority fermions is ex- 
pected. We note that our results have a direct relevance 
to the recently realized Fermi-Fermi mixture of and 
^Li because of their large mass ratio [8-10]. 

It is worthwhile to point it out that the effective in- 
teraction among heavy fermions immersed in the Fermi 
sea of light fermions is an analog of the Casimir force 
among objects placed in a vacuum [16]. The origin of 
the Casimir force can be traced back to the modification 
of the spectrum of zero point fluctuations of the electro- 
magnetic field. In our case, the role of vacuum is played 
by the Fermi sea of light fermions whose spectrum is 
modified by the presence of heavy fermions [17]. 

This paper is organized as follows. In Sec. II, we deter- 
mine the effective interaction between two heavy fermions 
induced by the interaction with the Fermi sea of light 
fermions as a function of the s-wave scattering length 
using the Born-Oppenheimer approximation. (The effec- 
tive interaction for a general number of heavy fermions 
is derived using the functional integral method and eval- 
uated in Appendix A.) With having in mind the appli- 
cation to the p-wave pairing of heavy fermions, we com- 
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pute the p-wavc projection of the effective interaction in 
Sec. III. Here we consider both cases where the heavy 
fermions are in three dimensions (3D) and where they 
are confined in two dimensions (2D) while hght fermions 
are always in 3D. The latter case corresponds to the two- 
species Fermi gas in 2D-3D mixed dimensions studied in 
Refs. [18, 19]. Finally summary and concluding remarks 
are given in Sec. IV. 



II. CASIMIR INTERACTION BETWEEN TWO 
HEAVY FERMIONS 

A. Born-Oppenheimer approximation 

Consider a two-species Fermi gas with unequal masses 
M and m interacting by an interspecies short-range inter- 
action. When the mass ratio is large M/m ^ 1. we can 
employ the Born-Oppenheimer approximation to com- 
pute the interaction energy of heavy fermions immersed 
in the Fermi sea of light fermions. Here we concentrate on 
the case of two heavy fermions. The formula for a general 
number of heavy fermions is derived in Appendix A using 
the functional integral method. The wave function of a 
light fcrmion interacting with two heavy fermions fixed 
at positions Xi and X2 satisfies the Schrodinger equation 
(here and below h= 1): 



E^{y;xi,X2) 



(1) 



The interspecies short-range interaction between the 
light and heavy fermions is taken into account by impos- 
ing the short-range boundary condition; ^{y — > Xi) oc 
l^j'"^ I — i + 0(|y — aJil), where a is the s-wave scattering 
length. Below we determine the energy spectrum of the 
light fermion both for bound states and continuum states 
as a function of a and the separation between the heavy 
fermions r = Xi — X2- 
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For the bound state E = —-^ < 0, the solution of 
Eq. (1) takes the form 



p-K\y-xi\ p-K\y-X2\ 

^±{y) = -. ^±-. r. 



\y - xi\ 



\y - X2\ 



(2) 



where ± is the parity of the wave function. The upper 
(lower) sign corresponds to the even (odd) parity. From 
the short-range boundary condition, we obtain the equa- 
tion 



1 

a 



(3) 



The solution to this equation exists when \r\/a > =f1 and 
is given by 



a r V / 



(4) 



where W{z) is the Lambert function satisfying z = 
W{z)e^^^\ Accordingly we find the following binding 
energy depending on a and |r|: 



binding 



2m 



(5) 



for \r\/a > =f1- 

For the continuum state E = > 0, the solution of 

Eq. (1) takes the form 



sin {k\y — Xi\ 



S±) ^ sin {k\y-X2\ + S±) 



*±(y) = 

(6) 

where < J± < tt is the phase shift. From the short- 
range boundary condition, we obtain the equation 



k cos 5± ± 



sin (A:|r| + 6-^ 



sin^a 



The solution to this equation is easily found to be 

k\r\ ± sin(A;|r|) 



tanJi = — ■ 



i^±cos(fc|r|) 



(7) 



(8) 



We then suppose that the system is confined in a large 
sphere with a radius i? ^ |r| and the two heavy fermions 
are located near its center; xi = ^ and X2 = By im- 
posing ^±(|y| — > -R) — > at the boundary of the sphere, 
the momentum k is discretized as 



k^R + 5+ = riTT, 
k-R + 6- = 



TT 



(9) 



with n = l,2,3, Therefore the energy level becomes 



2m 



(10) 



Now we fill the above-obtained energy levels with the 
light fermions. The total energy of such a system is given 

by 



E^ 



2m 



N 

E 

n=l 



2m 



(11) 



where 2A'' is the number of interacting light fermions in 

the continuum states. We compare the energy of the 
interacting system with that in the noninteracting limit. 



N 



E, 



free 



= E 



71=1 



2m 



(12) 



5±=0 



and define the energy reduction as AE'drj) = E — Eircc- 
Then we take the thermodynamic limit N, R oo with 
the Fermi momentum of the light fermions = Nw/R 
fixed. If we notice that the summation over n is dom- 
inated by n ~ ^ 1 and neglect small corrections of 
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FIG. 1; (Color online) Interaction energy of two heavy fermions v{kFr) as a function of their separation kpr. Three different 
curves correspond to the BCS regime (afcp)~^ = — 5 (dashed curve), unitarity limit (afcp)~^ = (solid curve), and BEC regime 
{akp)~^ — 5 (dotted curve). The right panel is a magnification of the left panel. 




FIG. 2: (Color online) Evolution of the interaction energy v(kFr) as a function of the inverse s-wave scattering length (afep)"^. 
Left panel: (afcp)"^ is varied from —2 (curve with the smallest amplitude) to (largest amplitude). The tiny oscillatory 
behavior existing in the BCS regime grows toward the unitarity limit and makes a small hump at r ~ fcp^. Right panel: 
(afep)~^ is varied from (curve with the smallest height) to 2 (largest height). The hump further grows in the BEC regime 
and develops a repulsive core at r ~ a. 



0{1/N), we obtain the following simple expression for 
the energy reduction: 



AE{\r\) = 



2m 







rriTT 



We note that the phase shifts given in Eq. (8) are defined 
to be in the range < S±{k) < tt. The same result 
can be obtained on a more general footing by using the 
functional integral method (see Appendix A). 



B. Effective interaction between two heavy 
fermions 



When the two heavy fermions are infinitely separated, 
the energy reduction approaches 



AE{\r\ ^ oo) ^ 2/Xsi„gie, 



(14) 



where /Xsingie < is the chemical potential of a sin- 
gle heavy fermion immersed in the Fermi sea of light 



fermions [20]: 



A^singlc 



kp ak-p + [l + (aA:F)^] [f -l- arctan(aA:F) ^] 



2m 



7r(aA:F)^ 



(15) 



The energy reduction with 2/isingio subtracted is regarded 
as the interaction energy of the two heavy fermions: 



y(|r|)^i?(|r|)-ii;(|r|^oo) 
= A£;(|r|)-2/isi„gic. 



(16) 



y(|r|) represents the effective interaction between two 
heavy fermions induced by the interaction with the Fermi 
sea of light fermions. We note that the chemical potential 

of light fermions /i; = is fixed here instead of their 
particle number. 

It is convenient to measure the interaction energy V{r) 
in units of the Fermi energy of light fermions and intro- 
duce a dimensionless function v[kYr) as 



V{r) = ^v{kpr). 
Zm 



(17) 
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v{k-pr) is a function of the separation between the heavy 
fermions kpr and also the s-wave scattering length afcp- 
v{kFr) for three typical values of (afcp)""^ = —5, 0, 5 are 
plotted in Fig. 1. One can see the smooth evolution of 
the effective interaction between the two heavy fermions 
as a function of (ak-p)^^ in Fig. 2. In the BCS regime 
(akp)^^ < —1, the effective interaction is attractive at 
r < |a| and has a tiny oscillatory behavior at r > \a\. 
This oscillatory behavior grows toward the unitarity limit 
(afcp)"^ ~ and makes a small hump at r ~ fcp^ (see 
the left panel of Fig. 2). This hump further grows in 
the BEG regime (afcp)""'^ ^ 1 and eventually develops a 
repulsive core at r ~ a (see the right panel of Fig. 2). 

It is worthwhile to compare our nonperturbative result 
(16) with the perturbative calculation of the effective in- 
teraction in the BCS limit a —0. In the BCS limit 
we have k± ^ and the phase shift 



|a| <r,fcF^ 



d±(k) — > —[kr± sin(A:r)] 
r 

± I cos(A;r) [kr ± sin(fcr)] + 0{\af). 



Thus we can find the interaction energy to be 



(18) 



v{kFr) 



^2 2fcFrcos(2fcFr)-sin(2fcFr) , ^,|„|3^ 

(19) 

which has the same form as the Ruderman-Kittel- 
Kasuya-Yosida interaction between magnetic impurities 
in a Fermi liquid [21]. Accordingly the Fourier trans- 
form of the effective interaction V^(|r|) in the BCS limit 
becomes 



/ 



V{\p\)= dre-^P-^V{\r\) 



(20) 



2m 



2 + 



2fc, 



-f-MVn 
p\ ^kp J 



1 + 


IpI 




2kF 




1 - 


IpI 




2/sF 





which correctly reproduces the effective interaction to the 
leading order in a obtained in Ref. [12] ^ as it should be. 
Our nonperturbative result (16) can go beyond the per- 
turbative BCS regime to the strongly interacting unitary 
regime in a controlled way by utilizing M/m 1. 

Analytic expressions of v{kFr) in various limits are ob- 
tained in Appendix B. 



III. P-WAVE PROJECTION OF THE 
EFFECTIVE INTERACTION 

A. Lessons from BCS theory 

The physics of heavy fermions immersed in the Fermi 
sea of light fermions is described by the Hamiltonian 



dxdy tljl{x)tljl{y)V{\x - y\)iph{y)i^h{x) 



(21) 



where /i/j is the chemical potential of heavy fermions mea- 
sured from /^single- Strictly speaking, the pairwise inter- 
action V(|r|) obtained in the previous section [Eq. (16)] 
is valid only if there are two heavy fermions in the system. 
This is because the Casimir interaction energy is not pair- 
wise additive and also because if heavy fermions have a fi- 
nite density, it will perturb the Fermi sea of light fermions 
and thus affect the effective interaction V^(|r|). However, 
in the dilute limit ^/2M fih <C fcF, we expect that the 
above Hamiltonian is a good approximation to the sys- 
tem of heavy fermions because the probability to find a 
third heavy fermion near the two heavy fermions becomes 
small in the dilute system [22]. Indeed we will confirm 
in Appendix A that the Casimir interaction among more 
than two heavy fermions at large separations can be eval- 
uated quite accurately as a sum of pairwise interactions 
between each of the two heavy fermions. This result also 
supports the use of our Hamiltonian (21) to describe the 
physics of dilute heavy fermions immersed in the Fermi 
sea of light fermions. 

Another issue to be addressed is the instability of the 
system. When the mass ratio M/m exceeds the critical 
value 13.6 in pure 3D [23] or 6.35 in the 2D-3D mix- 
ture [18], the system will not be stable due to the Efimov 
effect. Therefore we need to assume the mass ratio to be 
large but smaller than the above critical value. For such 
a mass ratio, the effective interaction l^(|r|) obtained in 
the Born-Oppenheimer approximation is no longer exact 
but can be considered as a good approximation to the 
exact one. 

One of predictions we can derive from the Hamiltonian 
(21) is a pairing between the heavy fermions. The pairing 
gap of heavy fermions Ap is defined to be 



(22) 



For unequal masses, the effective interaction in the BCS limit 

becomes y(|p|) = - [ WM)a j^ ^i(^), where L{z) = 

1 + "'"T^ InNJ-ti is the static Lindhard function. The limit of 

2 42 I 1 — 2 

large mass ratio M/m ^ 1 coincides with Eq. (20). 



where iph{p) is the Fourier transform of iph{x) and 
y(|p— ql) is the Fourier transform of y(|r I) [see Eq. (20)]. 
Because of the Fermi statistics of heavy fermions, the 
pairing gap has to have an odd parity; A_p = — Ap. 
The standard mean-field calculation leads to the self- 
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FIG. 3: (Color online) P-wave effective interaction between two heavy fermions in three dimensions as a function of (afcp)"^. 
■^NhVe^i{p) in d = 3 is plotted for p/kp = 0.4, 0.8, 1.2, and 1.6. The dashed curves are perturbative results at leading order 
valid in the BCS limit (afcp)^^ — cx). 



consistent gap equation 

A^^-Jj^V{\p-q\)^[l-2n^iE,)]. (23) 



Here Ep = y \^-^ — f^hj + l^pp is the quasiparticle 

energy and n-p{Ep) — 1/ [eP'p^^ + l) is the Fermi-Dirac 
distribution function at temperature T . 

When the coupling between different partial waves can 
be neglected, we can solve the gap equation (23) easily 
by using the weak-coupling approximation ^. For a given 
orbital angular momentum i (odd integer), the pairing 
gap and the critical temperature are given by [24] 

|Ap|~Te^M/«expf-i^y (24) 



^ The weak-coupling approximation even at unitarity is justified in 
the dilute limit; y/2MJi^ < fcp. Prom Eqs. (B2) and (B3), the 
interaction energy at the mean interparticle distance of heavy 

fermions r = {2Mfj.h)~ 2> fep becomes V{r) ~ 2m'prTp' ^'^'^ 

is parametrically small compared to the kinetic energy K{r) 
1 

2Mr^ • 



Here Nh is the density of states of heavy fermions at the 
Fermi surface and Vi is the partial-wave projection of the 
effective interaction V{\p— q\) with |p| — \q\ — \/2Mjih 
fixed on the Fermi surface of heavy fermions and assumed 
to be attractive < 0. Because the lowest partial wave 
for identical fermions is ^ = 1, the dominant pairing is 
expected to occur in the p-wave channel. With having 
in mind the application to the pairing of heavy fermions, 
we compute N^Vi as a function of the s-wave scattering 
length akp from Eq. (16). Below we consider both cases 
where the heavy fermions are in three dimensions (d = 3) 
and where they are confined in two dimensions {d = 2). 



B. Heavy fermions in d = 3 

When the heavy fermions are in three dimensions d = 
3, the partial-wave projection of the effective interaction 
V{\p — q\) with IpI = \q\ = p fixed is given by 

1 

Vt{p)=^^J dcosePi{cos9)V{\p^q\), (25) 

where cos 9 = p ■ q. Multiplying it by the density of 
states of heavy fermions Nh{p) = ^fi, we obtain the 
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FIG. 4: (Color online) P-wave effective interaction between two heavy fermions confined in two dimensions as a function of 



fNhVe^i{p) in d = 2 is plotted for p/kp = OA, 0.8, 1.2, and 1.6. The dashed curves are perturbative results at 



leading order valid in the BCS limit (afcp) ^ —>■ — oo. 



dimensionless function representing the effective interac- 
tion between two heavy fermions for the given partial 
wave i: 



NhVfip) 



d cos 6 Pi (cos ( 



M /fcp 
27rm \ p / j_i 

z sin(z\/2 — 2 cos( 



X I dz 
Jo 



V2 - 2cos^ 



fcpZ 



P 



(26) 



When the projected effective interaction with p = 
y/2MjIh being on the Fermi surface of heavy fermions 
is attractive NhVi < 0, the pairing of the heavy fermions 
is expected to occur with the pairing gap and the critical 
temperature given in Eq. (24). 

Figure 3 shows jjNhVe{p) in the p-wave channel £ = 1 
as a function of (ak-p)^^ for four values of p/kp — 

0. 4, 0.8, 1.2, 1.6. We note that what is plotted is not 
the p-wave effective interaction NhVi{p) cx itself but 
^NhVi{p) that is independent of the mass ratio M/m 

1 . We can see that the p-wave effective interaction is at- 
tractive in the BCS limit (akp)~^ ^ — 1 being consistent 
with the perturbative prediction [12]. However, it turns 
out that the perturbative result has only a small range 
of validity. Our result shows that the p-wave attraction 



is generally weaker than the extrapolation of the pertur- 
bative calculation and, in particular, the p-wave effective 
interaction becomes a strong repulsion near the unitarity 
limit. This may be understood because of the repulsive 
hump of v{kFr) developed near the unitarity limit (see 
Figs. 1 and 2). We also find that the p-wave attraction is 
stronger for the larger value of p/fcp because the contri- 



I ■ bution from the attractive part of v{k-pr) 



at short distance r ^ [see Eq. (Bl)] becomes more 
significant for larger p/k-p- The maximum attraction for 
each value of p/k-p = 0.4, 0.8, 1.2, 1.6 is achieved around 
[akp)^^ « —2.7, —2.5, —1.7,-1.0, respectively. At such 
a value of the s-wave scattering length, the maximal p- 
wave pairing of heavy minority fermions immersed in the 
Fermi sea of light fermions is possible while the pairing 
gap and the critical temperature would be very small be- 
cause of the numerically weak attraction NhVi^i{p) ^ 1. 
One can perform the same analysis for the /-wave chan- 
nel {(. — 3) but will find even weaker attraction. 



C. Heavy fermions in d = 2 

When the heavy fermions are confined in two dimen- 
sions d = 2, the partial-wave projection of the effective 
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interaction — q|) with \p\ = \q\ = p fixed is given by 



V,{p) = - / decosi£e)V{\p-q\), (27) 



where cos 9 = p ■ q. Multiplying it by the density of 
states of heavy fermions Nh{p) = ^) we obtain the di- 
mensionless function representing the effective interac- 
tion between two heavy fermions for the given partial 
wave i: 



NhVeip) = 



M fkp 



2-Km \ p 

X 



Jo 



de cos(£6i) 



J dzz Jo(^z\/2- 2 cosOj ^(^^^ 



(28) 



In the case of c? = 2, the scattering length a should be 
regarded as the effective scattering length Ugs introduced 
in Ref. [18]. 

Fi gure 4 shows N^V^ij)) in the p-wave channel £ = 1 

as a function of (ak-p)"^ for four values of p/kp = 
0.4, 0.8, 1.2, 1.6. One can sec the similar behavior to 
the case of d = 3. Again the p-wave effective interaction 
is attractive in the BCS limit {akp)~^ <C —1 being con- 
sistent with the perturbative prediction [19]. Compared 
to the case of = 3, we find that the perturbative result 
has a wider range of validity. In the case of d = 2, the p- 
wave attraction is found to be stronger than that in = 3 
for the same value of p/fcp and the p-wavc effective inter- 
action becomes only a weak repulsion near the unitarity 
limit. This will be understandable because the contri- 
bution from the attractive part of v{k-pr) —c^/{k-pr)'^ 
at short distance r ^ [sec Eq. (Bl)] is less suppressed 
by the phase space factor and thus more significant in 
d = 2. Our result also shows that the p-wave attraction 
is stronger for the larger value of p/kp and the maxi- 
mum attraction for each value oip/kp = 0.4, 0.8, 1.2, 1.6 
is achieved around (afcp)"^ ~ —1-2, —1-2, —0.85, —0.52, 
respectively. At such a value of the s-wave scatter- 
ing length, the maximal p-wavc pairing of heavy minor- 
ity fermions in two dimensions immersed in the three- 
dimensional Fermi sea of light fermions is possible and 
the pairing gap and the critical temperature will be larger 
than those in the case of d = 3. 



the p-wavc projection of the effective interaction is attrac- 
tive in the BCS limit being consistent with the pertur- 
bative prediction, while it turns out to be repulsive near 
the unitarity limit. We found that the p-wave attraction 
reaches its maximum between the BCS and unitarity lim- 
its, where the maximal but still weak p-wave pairing of 
heavy minority fermions is possible. 

We also investigated the case where the heavy fermions 
are confined in two dimensions corresponding to the two- 
species Fermi gas in the 2D-3D mixed dimensions [18, 
19]. Because the p-wave attraction between the heavy 
fermions in two dimensions was found to be stronger 
than that in three dimensions, we expect the larger p- 
wave pairing gap of heavy minority fermions. Such a 
p-wavc pairing in two dimensions is especially interesting 
because the resulting system has a potential application 
to topological quantum computation using vortices with 
non-Abelian statistics [25, 26]. 

Although our controlled analysis is restricted to the di- 
lute heavy fermions in the limit of large mass ratio, our 
results indeed shed light on the phase diagram of asym- 
metric Fermi gases with unequal densities and masses 
(and spatial dimensions), in particular, in the strongly 
interacting unitary regime. Our results also have a direct 
relevance to the recently realized Fermi-Fermi mixture of 
''"K and ^Li because of their large mass ratio [8 10]. It 
will be an important future problem to study corrections 
to our results when the restrictions of the large mass ratio 
and the dilute limit of heavy fermions are relaxed. 
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APPENDIX A: CASIMIR INTERACTION 
AMONG HEAVY FERMIONS FROM 
FUNCTIONAL INTEGRAL METHOD 

1. Derivation of Casimir interaction 



IV. SUMMARY AND CONCLUDING 
REMARKS 

We have investigated a two-species Fermi gas with 
a large mass ratio interacting by an interspecies short- 
range interaction. Using the Born-Oppenheimer ap- 
proximation, we determined the interaction energy of 
two heavy fermions immersed in the Fermi sea of light 
fermions, which is an analog of the Casimir force, as a 
function of the s-wave scattering length. We showed that 



Here we derive the Casimir interaction among heavy 
fermions immersed in the Fermi sea of light fermions us- 
ing the functional integral method [27]. The action that 
describes the light fermions with the chemical potential 
fil interacting with the heavy fermions fixed at positions 
Xi is 



(Al) 



5 = y drdxipjir, x) ^'■('^^ 

~9oY^ jdTll>}{T,Xi)'4)l{T,Xi), 
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where r is an imaginary time. The partition function 
Z = e~^^ is given by 

Z = j Vil)iDil)\e-^ . (A2) 

We first insert the following identity into the integrand; 

1 = JJ / 1>riiDfii 5[i}i{T,Xi) - ■qi{T)]5[i)\{T,Xi) - f?j(r)] 

i 

(A3) 

and then exponentiate the delta functions by introducing 
auxiliary fields, 

1 = n / 'Dvi'DfjiVatVai (A4) 

Now the partition function can be written as 

Z= I Vi^iVij\^Vf^,Vfi{Da^Va^e-^' , (A5) 

where the action becomes 

S' = J drdx ^1)\{t, x) ~ ^ ~ /^'^ ^li^^ ^) 

-9oYl J dr f)i{T)rii{T) 

_ (A6) 
-izL J dTai{T)[ipi{T,Xi) - r/i(r)] 

i 

I dT[il}\{T,Xi) - fii{T)]ai{T). 

i ■' 

We can easily integrate out ipi and ip\ fields and rj and 
f) fields to lead to the partition function 



Z = Zq j l[VaiVai e"^", 



(A7) 



where Zq = e~^^'""' is a partition function of noninter- 
acting light fermions and the action S" in the momentum 
space becomes 



^„ ^ f dwdp _ e^P^-,) 



-aj{Lv) 



(AS) 



Finally, by integrating out a and a fields, we obtain the 
following expression for the partition function: 



Z = Zq exp 



/3 [ ^\ndet Mi j{iuj) 

J 277 



(A9) 



where My is the scattering matrix given by 



90 



dp 



= < 



1 

m 



(AlO) 

\/—2miuj — 2mfii^ for i = j 



, 27r|a;i — a;^ 

Here wc introduced the s-wave scattering length a 
through 



1 r dp 2r 
go J (27r)3 p 



2m 



m 

2-Ka' 



(All) 



Therefore the reduction in the grand potential O com- 
pared to that in the noninteracting limit Q-hee is 

I^n{{xi}) = ~ I ^\ndetMij{iuj). (A12) 

J 2lT 

By subtracting Af2 in which all heavy fermions are in- 
finitely separated, the interaction energy of the heavy 
fermions becomes 



duj ^ det Mij (iuj) 
dct Mf^{iuj)' 



(A13) 



where M!j^ is a diagonal scattering matrix composed of 
Mjj. This is the generalization of V(|r|) in Eq. (16) to a 
general number of heavy fermions fixed at positions cc,. 



2. Casimir interaction between two heavy fermions 

It is worthwhile to reproduce the result obtained in 
Sec. II in the case of two heavy fermions. By deforming 
the path of the integration over iu = iiv + m into 



^ + / 

-ioo-\-fj,i J — oo— iO+ J fj, 



(J-l p — oo-\-iO'^ 

+ 



(A14) 



we obtain the following expression for the grand potential 
reduction: 



diy 



1 



A^l{\r\) = i / — In - - V-2mv + 
277 \ a 



- yj —2v(iv\r\ 



+ ^y-ln(^--./=2^ 

(A15) 

with r — Xi — X2- Wc separate the integral into the con- 
tributions from bound states Re[i^] < and continuum 
states Re[i/] > 0. For the bound state contribution, we 
pick up the singularity in the integrand (A15) as 



dv 



1 



I I T^^^l \/—2mv ± 

/ReM<o 27r \ a 



-^J — 1vav\r\ 



= i/± + const. 



(A16) 
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FIG. 5: (Color online) Interaction energy of three heavy fermions fixed with the same separations r. vs{kFr) as a function of 
fcpr is plotted for (afcp)~^ = — 1 (left), (middle), and 1 (right). The dashed curves represent the sum of pairwise interaction 
energies 3v2{kFr). 



V4(kfr) 

O.lOp 




FIG. 6: (Color online) Interaction energy of four heavy fermions fixed with the same separations r. ii4(fcpr) as a function of 
fcpr is plotted for (afcp)~^ = —1 (left), (middle), and 1 (right). The dashed curves represent the sum of pairwise interaction 
energies 6w2(fcFr). 



where < is the binding energy satisfying 

I g-\/-2mi^|7-| 

- - y/-2miy± ± =^ 0. 



(A17) 



For the continuum state contribution, defining the phase 
shift by 



n — 6± (i^) 

= arg I - — ^ —2mv — iO+ ± 



(A18) 



-V-2mjy-i0+|r| ' 



the integral in Eq. (A15) can be written as 

■/ 



dv I I 



II — In \/ —2mv ± 

[ii- \ dv ■ 



— \/—2mv\r\ 







TT 



(A19) 



Therefore, after dropping the unimportant constants, we 
find that the grand potential reduction in the case of two 
heavy fermions is given by 



An{\r\) = iy+ + u_ - dv 
Jo 



d+iv)+S^iv) 



(A20) 



This result is equivalent to Ai?(|r|) in Eq. (13) with v± = 

~5m' ~ ^^'^ ~ Im thus provides the same 
interaction energy T^(|r|) as Eq. (16). 



Casimir interaction among three and four heavy 
fermions 



We now evaluate the Casimir interaction energy (A13) 
among three and four heavy fermions fixed with the same 
separations \xi — Xj\ = r. We measure the interaction 
energy in units of the Fermi energy of light fermions 

VN{r) = ^VNikpr) with the subscript TV indicating the 
number of heavy fermions. The dimensionless functions 
VNikpr) with three typical values of (afcp)"^ = 0, 1 
are plotted as functions of kpr in Fig. 5 for = 3 and 
in Fig. 6 for A^ = 4, together with the sum of pairwise 
interaction energies — —V2{kFi') for comparison. 



We can see that the Casimir interaction among more 
than two heavy fermions can be reproduced quite accu- 
rately by the sum of pairwise interactions between each of 
the two heavy fermions, in particular, at long distances, 
although the behaviors at short distances are slightly 
overestimated. A similar observation has been made in 
Ref. [17] in the system of hardcore spheres immersed in a 
background Fermi sea. These results support the use of 
the Hamiltonian (21) only with the pairwise interaction 
to describe the physics of dilute heavy fermions immersed 
in the Fermi sea of light fermions. 
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APPENDIX B: INTERACTION ENERGY IN 
VARIOUS LIMITS 



In particular, in the unitarity limit akp —^oo,we find 



Here we evaluate the effective interaction between two 
heavy fermions v{k-pr) in Eq. (17) in various Hmits where 
analytic expressions are available. 

At short distance r <C a, fcp ^, we obtain 



v{k-pr) 



2c 1 



1 + c + 1 



+ 



{kFrf l + cak^r {I + cf a^k§ 
2akF + [l + {akp)^] [n + 2 arctan(afcF)"^] 



7r(afcF)^ 



-l + 0(r), 



(Bl) 



where c = 0.567143 is a solution to c = e^'^. 

On the other hand, at long distance r a, fcj^, we 
obtain 

vikpr) (B2) 

2(aA;F)3sin(2fcFr) - {akpf [{akpf - l] cos(2fcFr) 
^ 7r[(aA:F)2 + l]'(fcFr)3 

4:{akF f [{akp)^ - l] cos(2fcFr) 
27r[(afcF)2 + 1]^ (fepr)^ 
_ (gfcF)- [(ah-)^- - + l] siii(2/vi./0 

27r [(afcF)2 + 1]^ (A;Fr)4 



+ 0(r-5). 



v{kFr) — > — 



+ 



cos(2fcF'') sin(2A;Fr') 
7r(fcFr)3 ~ 27r(A;Fr-)'* 
2 cos(2fcFr) + cos(4fcFr) 
47r(fcFr")^ 



(B3) 



0[{kFry 



In the BCS or BEG limit \a\ <C r, fcp ^, we obtain 



v{kFr) (akp) 



22kFr cos{2kFr) — sin(2fcFr) 
27r(fcFr)4 



+ (c/cf) 
+ 0{a^). 



3 2A:f;-c-{),s(2Av/-) + [2(/vf;-)- - l] riiu(2fcFr) 



7r(fcFr)^ 



(B4) 
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